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1 Introduction 

In a recent paper Q|, we have started the NLO calcula- 
tion of the jet vertex which represents one of the build- 
ing blocks in the production of Mueller-Navelet jets || 
at hadron hadron colliders and of forward jets || in deep 
inelastic electron proton scattering. Both jet production 
processes are of particular interest for studying QCD in 
the Regge limit or in the small- x limit: they provide kinc- 
matical environments for which the BFKL Pomeron pj 
applies. Previous experience shows that existing leading 
order calculations |||||7|] are not accurate enough to allow 
for a reliable comparison with experimental data. NLO 
calculations are available for the BFKL Pomeron 
but consistent next-to-leading order analysis of data at 
the Tevatron, at LHC, or at HERA require the NLO cal- 
culations also of the jet vertex and of the photon impact 
factor. As to the jet vertex, in our previous paper j[| we 
have presented the first part, namely the quark-initiated 
jet vertex. In the present paper, we complete the NLO 
analysis of the jet vertex with the gluon-initiated part. 
The NLO calculation of the photon impact factor is being 
pursued by two independent groups [jl^,|l). As a result of 
these combined efforts, it will be possible perform consis- 
tent NLO studies of BFKL predictions in hadron hadron 
colliders, in ep deep inelastic scattering, and in e + e~ scat- 
tering processes. 

A particular theoretical challenge in computing the jet 
vertex for the processes mentioned before is related to the 
special kinematics. The processes to be analyzed is illus- 
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trated in Fig.|j]: the lower gluon emitted from the hadron 
H scatters with the upper parton q and produces the jet J. 
Because of the large transverse momentum of the jet, the 
gluon is hard and obeys the collinear factorization. In par- 
ticular, its scale dependence is described by the DGLAP 
evolution equations fl2| . Above the jet, on the other hand, 
we require a large rapidity gap between the jet and the 
outgoing parton q: this kinematic requirement is described 
by BFKL dynamics. Consequently, the jet vertex lies at 
the interface between DGLAP and BFKL dynamics. As an 
essential result of our analysis we find that it is possible 
to separate, inside the jet vertex, the collinear infrared 
divergences that go into the parton evolution of the in- 
coming gluon from the high energy gluon radiation inside 
the rapidity gap which belongs to first rung of the LO 
BFKL ladder. In @ this was demonstrated for the quark- 
initiated vertex, and in the present paper we present the 
generalization to the gluon part. 

Our paper will be organized as follows. In order to 
make the reading as convenient as possible, we first re- 
view the general framework, which will be the same as 
in our previous paper Sections || and || then contain 
the virtual and real corrections, resp., and in section [5] we 
combine both results. Section [| contains, as a consistency 
check, a comparison of our jet vertex with the gluon im- 
pact factor. In section we use our results to define the 
general NLO jet production cross section. The concluding 
section contains a short summary and an outlook of future 
steps. 



2 High energy factorization 
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2.1 General framework 

We describe the kinematics of the hadron (H) quark (q) 
collision in terms of light cone coordinates 



P 1 * = (P + ,P ,P) , P 



± _ p°±p 3 



(i) 



where the light-like vectors pu and p q form the basis of 
the longitudinal plane: 



-,0,0) , s:=(p H + Pq y (2a) 



(2b) 
(2c) 



Ph 



P,= (0,\/2>° 



e Vi e -Vi 

Pi = Ei[ —?=rAi 



In the last equation we have introduced a parameteriza- 
tion for the i-th particle in the final state in terms of the 
rapidity yi (in the ph + Pq center of mass frame), of the 
transverse energy E% = |p 4 1 and of the azimuthal unit vec- 
tor fa \\ Pi . 




t 

Fig. 1. Diagrammatic representation of the high energy pro- 
cess studied in this paper. H is the incoming hadron providing 
a gluon g with distribution density /; q is the other incoming 
particle which will be taken to be a quark; J denotes the jet 
produced in the forward direction (w.r.t H) and i is the generic 
label for outgoing particles. 



According to the parton model, we assume the physi- 
cal cross section to be given by the corresponding partonic 
cross section dcr (computable in perturbation theory) con- 
voluted with the parton distribution densities (PDF) / a of 
the partons a inside the hadron H . A jet distribution Sj 
selects the final states contributing to the one jet inclusive 
cross section that we are considering. 

In terms of the jet variables — rapidity, transverse 
energy and azimuthal angle — the one jet inclusive cross 
section initiated by the gluons in hadron H can be written 
as 



da 
dJ 



da qH 



dyjdE, 



= I dxda qg (x)Sj(x)fW(x) . (3) 



The incoming gluon carries a fraction x of the longitudinal 
momentum of H, while its transverse motion is negligible 
in the high energy regime: 



Pg = xp H 



■rd -,0.0 



(4) 



In our analysis we study the partonic subprocess g + q 
X + jet in the high energy limit 



A 



QCD 



< E 2 j - -t (fixed) < s -> oo 



(5) 



2.2 The Jet vertex at lowest order 

In the high energy regime (^|), the lowest order (LO) con- 
tribution to the jet cross section is dominated by gluon 
exchange in the i-channel, as shown in Fig. FJ. 

The corresponding expression has been already ob- 
tained in and is given by 



= dx dk h q °Hk)V g W (k,x)ff\x) (6) 



da(°) 
dJ 

in terms of the LO jet vertex 

V g (°)(fc,x):=/4°)(fc)4 2) (k;x) , (7) 
which is just the product of the LO gluon impact factor 

2 1+£ a 



hf){k):=M%, JV = — 



(8) 



and of the jet distribution for two particles in the final 
state 

Sf\k;x) := sf\pi,p 2 ;p g ,p q ) 



Xj : 



Ejey j 



(9) 



2.3 Ansatz for the factorization formula 

According to the analysis of Ref. [Q, we propose a high 
energy factorization formula for the description of the high 





7 > 







Fig. 2. Leading diagram at lowest order for gluon-quark scat- 
tering: the interaction occurs via gluon exchange in the t- 
channel. 
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energy quark-hadron interaction with a jet in the final 
state. In this formula a quark impact factor h q , the gluon 
Green's function G, the gluon-initiated jet vertex V s and 
the gluon PDF f g are convoluted in both transverse and 
longitudinal variables to provide the jet differential cross 
section as follows: 



da 



dJ 



= / dx 



dk dk' h q (k)G(xs, k, k )V g (k , x)f g (x) 

(10a) 

G(xs 7 k. k') := S(k - k') + a s K^\k, k') log — + 0(a 2 s ) . 

(10b) 



•so 



The gluon Green's function contains by definition all the 
energy dependence of the process and is given in terms of 
the BFKL kernel Q K (k, k'). By writing the perturbative 
expansions for the quark impact factor, the jet vertex, and 
the PDF 



h = + a s h^ - 
V = V<® + a s V^ 

/ = / (0) +«s/ (1) - 



(11a) 
(lib) 
(11c) 



our ansatz corresponds to the following structure for the 
one-loop cross section: 



da(D 
dJ 



dx / dk 



dk' 



h^(k)K^(k,k') 



xlog-V^(k\x)f^(x) 

+ h^(k)V^(k,x)f^(x) 
+ h^(k)V^(k,x)f^(x) 

+ h^(k)V g il Hk,x)ff\x) 



(12) 



which is obtained simply by expanding Eq. (10a) up to 
relative order a s . 

For the first order correction to the partonic impact 
factor, foW , which appears on the second line, we can use 
the known expression of Rcf . (ll|[l4| , and for the correction 
to the PDF, f( x \ we have the usual convolution with the 
LO Altarelli-Parisi splitting functions: 



*Jf(«-*)-£i(#)'?/.T P - ({) ^(f 



(13) 



Because of the definition ( |14[) for a s = g pb /At: 



log 47r)], Eq. (|13|) defines the one- loop PDF in the MS 
scheme. Finally, the correction term is what we want 
to compute in this paper. 

Equations ([H]) and ([ll]) constitute a highly non trivial 
ansatz, which will be shown to depend upon a careful sep- 
aration of singular and finite pieces. Our main task will 



consist to identify the collinear singularities ( |l3| ) to be ab- 
sorbed in the parton densities, to check cancellation of the 
remaining infrared singularities, and, finally, to separate 
the terms proportional to log s which belong into the first 
line of (|T2]). The remaining finite (in e) and constant (in s) 
term will eventually be interpreted as one-loop correction 

to the jet vertex, 



(i) 



3 Virtual corrections 

For the one-loop analysis of the gluon-initiated jet pro- 
duction process we adopt dimensional regularization in 
D = 4 + 2e dimensions and define, according to the MS 
scheme, the bare dimensionless coupling a s as a function 
of the dimensionful bare coupling g and of the renormal- 
ization scale /i as follows: 



= ct® := 



ffi 2e r{i 



(47T) 



l+e 



(14) 



The one-loop analysis of the virtual corrections can be 
carried out in the same way of the quark-initiated case 
presented in [ l) . Discarding all terms suppressed by powers 
oop quark-gluon cross section can be derived 
| and the ensuing virtual contribution to the 



of s, the one- 
from Ref. 
jet cross section reads 



= a 8 / dx J dk h«>\k) 



2w (1) (fc)log 



n q (k) + n g (k) 



xV^(k,x)6°Hx), 



(15) 



h q °\k) — Cp/CAhg (k) being the LO quark impact fac- 
tor. The first term represents the leading log s (LL) contri- 
bution to the virtual corrections. The coefficient of logs, 
namely 2a/ 1 ), constitutes the virtual part of the leading 
BFKL kernel and is just twice the one-loop Regge-gluon 
trajectory 



;«(fc) 



C A 1 r 2 (l + e) (k 



7T 2er(l + 2e) \fi' 



(16) 



The e-pole reflects the presence of a soft singularity which 
will be compensated by an opposite one in the real part 
of the kernel. 

The non logarithmic terms in Eq. ( p"5| ) represent the 
next-to-leading logs (NLL) contribution to the virtual 
corrections after the subtraction of the UV e-pole occur- 
ring in the renormalization of the coupling 



a s (n 2 ) : 



,(0) 



1 



,(o)! 



(17) 



They are expressed in terms of the virtual corrections to 
the impact factors denoted by 77. The virtual correction 
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to the gluon impact factor reads 



n g (k) 



1 


11 


5tt 2 


67\ 


Ca 


-— ■+ 


D£ 


~L2~ 


~ 36 y 


7T 


1 


5 N 




- b log 


fc 2 " 


~3e 




7T 


>_ 



(18) 



where &o = (HCU — 2Nf)/12ir is the first coefficient of the 
/3-function. Any occurrence of a s in Eq. ( |l5| ) and in all 
other coming formulae is to be understood as a s (/i 2 ). 

The gluon impact factor virtual correction ( Jig ) shows 
double and single poles in e. These poles are of IR origin 
and are due to both soft and collinear singularities. Partly 
they will cancel against the corresponding singularities of 
the real emission corrections, leaving a simple pole that 
will be absorbed in the redefinition of the PDFs. This will 
be shown in Sec. o. 



into two semispaces defined by y{, y' 3 > (lower half) and 
y' 3 < < y{ (upper half) Q. 

The form of the partonic differential cross section turns 
out to be quite simple when restricted to one of the two 
halves of the phase space. For the "lower half region" 
y[,y' 3 > 0, which corresponds to z > z cut '■ — the 
cross section can be cast into the general form 



(20) 

dS-qg-yft, - h q °\k)F fin (k,k',z)h^(k') dkdk'dz , 
(z > Zcut) , (21) 



where the function F depends on the particular final state. 
For quark-gluon scattering, we have two contributions. 
The qqq final state term is given by 



4 Real corrections 

The real emission corrections of the quark-gluon collision 
that we are considering, involves three partons in the final 
states: an outgoing quark (labelled by "2") — which is 
nothing but the scattered incoming quark q — and two 
additional partons ("1" and "3"), which can be either a 
qq pair or two gluons (see Fig. ||). 




Fig. 3. Labelling the two-to-three parton scattering process. 
The two dashed lines denote a qq pair or a gluon pair. Note 
the definition of the transferred momenta fc and fc'. 



F^(k,k',z) = ^-N f T J 



V qg (z,s) 1 



V qg (z,e) = l 



C F qk' 

2z(\- z) 
l + e ' 



(22) 



(23) 



tt e = ^T(l - e)^ , 



where q = fc — k' , Tr = 1/2 and Vg^z, e) is — apart from 
a missing Tr factor — the real part of the g — > q splitting 
function in 4 + 2e dimensions. 



The qgg contribution can be written as 



F qgg (fc,fc',z) = -^C A 



V m (z) 



V m (z) 



2ir 7T L 

z 2 k' 2 + (1 - z) 2 q 2 - z(l - z)q ■ k' 
q 2 (q — zk) 2 
l + z 4 + (l-z) 4 



(24) 



2z(\-z) 



(25) 



We parametrize the kinematics of the process in terms 
of Sudakov variables of two exchanged momenta: 



k := p q - p 2 = -wp g + wp q + , fcj. = (0, 0, fc) 



(19a) 

fc' := Pi — p g = -zp g + zp q + k' ± , k' ± = (0, 0, fc') . 

(19b) 

Note that the transverse energies introduced in Eq. ( |2c| ) 
correspond to E± = \k'\,E 2 = (fc),^ = \k — fc'|. 

The partonic differential cross sections has been com- 
puted (ll|[l4| for both cases in the high energy regime, 
where we neglect terms suppressed by powers of s. In- 
troducing the rapidity y' = y+ \ log - "measured" in the 
partonic center of mass frame, we can split the phase space 



where T , gg (z) is — apart from the missing colour factor 
2Ca — the real part of the g — * g splitting function (in 
any dimension). 

The real correction to the upper quark impact factor 
receives contribution only from the qgg final state in the 
"upper half region" of the phase space and from the vir- 
tual correction contribution of Eq. ( |l5| ) coming from the 
q impact factor correction 77 q . It can be computed in a 
similar way as presented in [Q and it will not be discussed 
further here. We simply report the final result adapted to 



1 By momentum conservation there is always a particle with 
y > 0, which can be labelled by "1" without loss of generality. 
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the case of incoming gluon: 



d(7^ <0) da^ 



dJ 



dJ 



x log 



dx J dfcdfc' ^(fe)if(°' real )(fe,fe') 

|^yf(fc',,)/f)(x) (26a) 

max(£/ 2 , £-3) 

dx f dk h^(k)V^ 0) (k,x)f^(x) , (26b) 



The structure of the infrared (IR) singularities for the 
gluon-initiated vertex, is somewhat more entangled than 
that of the quark-initiated one. This happens in particu- 
lar when the gluon pair is emitted in the final state. Since 
gluons are indistinguishable, each of the gluons can con- 
tribute to the LL term, and the LL subtraction needed 
to extract the jet vertex cannot be defined on the same 
footing as in the quark case 0. 

In the next sections we develop an explicit procedure 
for separating the LL contribution from the collinear sin- 
gular terms. We check the cancellation of the soft singular- 
ities and isolate the finite expression for the gluon-initiated 
jet vertex. 



4.1 Jet definition 

We begin with a brief review of the jet definition. Follow- 
ing the arguments given in jl6j, we require that the jet 

distribution Sj 1 ^ , selecting from a generic n-particle final 
state the configurations contributing to our one jet inclu- 
sive observable, be IR safe. This corresponds to the fact 
that emission of a soft particle cannot be distinguished 
from the analogous state without soft emission. Further- 
more, collinear emissions of partons cannot be distingui- 
shed from the corresponding state where the collinear par- 
tons are replaced by a single parton carrying the sum of 
their quantum numbers. 

Listing the possible IR singular configurations, only 
the emission of a gluon (1 or 3) with vanishing momen- 
tum gives rise to soft singularities. Collinear singularities 
arise in collinear emissions of partons that couple directly 
to each other. The list of all possible collinear singular 
configurations, in an obvious notation, reads as follows: 



enters the jet, and only particles 1 and 3 play a role in 
building up the jet. 

In terms of the variables p 1 , p 3 , p^ , p^ , the 3-particle 
jet distribution 



4 3) (pi, p 3 , = 4 3) k - k '> *) 

Ph Ph 



'■= Sf\pi,P2,P3;Pg,Pq) 



(28) 



has to satisfy the following properties in order to be IR 
safe (cfr. §): 



1 soft : 


S { f\o,p,x;x) 


= Sf\p;x) 


(29a) 


3 soft : 


sf } {p,O,0;x) 


= sf{p-x) 


(29b) 


1|| 3: 


Sf\(l-\)p,\p 7 \x;x) 


= Sf\p;x) 


(29c) 


gill: 


4 3) (0,p,C;x) 


= sf\p;0 


(29d) 


g||3: 


4 3) (p,0,t;x) 


= sf{p-x- 


0- 

(29e) 



Note that, since quark 2 does not participate in the jet, 
the collinear properties of the jet di strib ution are applied 
only to the configurations listed in (27a). 

In the following sections these relations will be used 
when extracting the divergences of the real emission. 



4.2 qq pair in the final state 

We have already observed in the previous section that 
there are two kind of final states, differing in the type 
of the partons denoted by 1 and 3: either a qq pair or a 
gluon pair. Here we consider the contribution of the final 
state with qq pair. In this case, the main contribution 
occurs when the qq pair emitted with a small invariant 
mass in the fragmentation region of the incoming gluon. In 
practice, we need only to consider the "lower half" phase 
space 2/3 < contribution. The corresponding Feynman 
diagrams are shown in Fig. ||. 

The starting formula is derived from Eq. (^), using 
Eq. ( |2"8| ) for the jet distribution and Eqs. (|2l]) and ( p2[ ) for 
the partonic cross section: 



do- 
ck/ 



qqq 



2tt 



N / dfcdfc' /4 0) (fc) 



2 II 3 



(27a) 
(27b) 



dz N f T R - 



V qg (z,e) 1 



fe' 2 t 



It is important to note that, in the kinematic regime we 
are considering, configurations in which quark 2 is emitted 
outside the fragmentation region of quark q are strongly 
suppressed, i.e., quark 2 never belongs to the jet produced 
in the forward direction of gluon g. This is because the 
propagator of the exchanged particle provides a suppres- 
sion factor ~ t/s with respect to the situation where quark 
2 is in the fragmentation region of quark q. Therefore, 
we can safely neglect the configurations in which quark 2 



Cf + Ca 



z(l — z)k ■ q 



(q-zk) 2 
J dx Sf\k',q,xz;x)f^(x) 



(30) 



Since the integrand in the RHS of Eq. (|3(j) is regular in the 
z — > limit, we can replace the lower limit of integration 
in z: z cut — > 0. This amounts to a negligible error in the 
high energy limit of order t/s. 
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Fig. 4. Feynman diagrams contributing to quark-gluon scattering in which the gluon fragments in a qq pair. 



Cf term 



Tr[z 2 + (1 — z) 2 ] one can write 



Let us first consider the Cf piece. In this term the two 
IR singularities of collinear origin (3 || g <^=^ q = 
or 1 || g <^=> fe' = 0) can be disentangled employing a 
simple fract decomposition: 



1 



1 



q 2 k' 2 
and we can write 



da Cp _ da^l 



dJ 



d.J 



fe' 2 



dJ 



1 



1 

fe 72 " 



/ a, „<»><*> j^mn) 



dJ J q J ir £ (k- 

where the integrand I is given by 



(31) 



(32) 



da™ 11 
dJ 



T = ^ N f J dk /l q° ) ( fc ) / dx 



1 fA 2 ^ - 

2 



, , , dzP qg {z)V^{k,xz)f^\x) 
£ VA 1 / Jo 

dzz(l-z)V<°\k,xz)fW(xj\. (36) 



The finite contribution of the Cp term can be computed 
at e = and reads 



dcr finite 

dJ 



^N f J dk h^(k) J^-J dxj dzP qg (z) 



h^(k' 



——^sf(k',q,xz;x) 



j^G(A 2 -k' 2 )VW(k,xz)\fW(x). (37) 



1 



2tt 



dz P q g(z, e) 



k' 2 + (fe - fe') 2 io 
/" da;^ 3) (fc',fc-fc',a;z;a;)/^(a;) . (33) 



(3) 

Because of the symmetry of V qg and S j under the q <-> q 
exchange, i.e., {z <-> 1 — z, fe' <-> fe — fe'}, it holds 7(fc'; fe) = 
/(fe — fe'; fe) and thus the two terms in Eq. ( |3^ ) are equal. 

The separation of the collinear singularity is performed 
by means of the subtraction method by decomposing 



Ca term 



In the Ca term there is only one final state collinear sin- 
gularity when 3 || 1, corresponding to q — zk — 0. Also in 
this case we separate the collinear singularity 



^c A da-; 1 da*f ° 



dJ 



dJ 



dJ 



(38) 



d^c* 
dJ 



dcr?? 11 
Cf 

dJ 



dfJ finite 

dJ 



(34) 



where the collinear singular term depends on an UV cut- 
off A and is defined by 



d<7 co11 
dJ 



coil f r Ah' 

C F := 2 J dk (fe) y ^1(0; k)0(A 2 fe' 



(35) 



In both collinear limits 1 || g and 3 || g the jet distribu- 



tion reduces to -► 5^ ; (fc; x(l - z)) (cfr. Eq. (|29dJ,e)) 



:< 2 



Introducing the gluon-to-quark splitting function P qg (z) 



by means of the subtraction method: we define the colli- 
near singular term as the residue of the l/(q — zk) 2 pole, 
integrated in a circle of radius A centered in the singular- 
ity. For simplicity we use the same cut-off A introduced 
in the Cf term. In this case, the jet distribution can be 



(3) 11 

simplified thanks to Eq. (|29c|) : Sj — > Sy j (k;x) and the 
z integration can straightforwardly be performed, yielding 



dJ 



dfc/4 0) (fc) / dxV R i0 \k,x)6°\x) 



1 



1 

12 



(39) 
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and 



^.finite 
dJ 



^-N f Jdkh^(k)JdxfP(x) 

x ( dzAs(z) /v(^w 



q 2 k 



,0(A 2 -{q- zk) 2 )sf ] {k;x) (40) 



4.3 Real corrections to the jet vertex from qgg final 
state 

The starting formula can be derived from Eq. using 
Eq. HI for the jet distribution and Eqs. (|l|) and @) for 
the partonic cross section. A more convenient expression 
can be written if one notes that the splitting function for 
gluons can be decomposed as 



T ge (z) = P(z) + P(l - z), P(z) = I - + - 1 (1 - z). 

(41) 

Using the symmetry of the partonic cross section and of 
the jet distribution under the exchange of the two gluons: 
{z <-> 1 — z,k' «-> k — k'}, the jet cross section can be 
rewritten as 

= 2i / dfc h®(k') / — h^(k) / dz P(z) 

dJ 7T J J TT £ J z cut 



2 k' 2 + (1 - z)q ■ (g - zk) 
q 2 (q — zk) 2 

x / dxSf\k',q,xz;x)f^\x) 



(42) 



It is important to observe that now the "splitting func- 
tion" P has a pole only at z = 0, while it is regular (ac- 
tually vanishes) at z = 1. The reason is that now we are 
employing an asymmetric treatment which causes gluon 3 
to be the only responsible of soft singularities and of the 
central region LL contribution. 

The expression ([42]) looks now pretty similar to the 
real contribution of the quark-initiated case. We can there- 
fore proceed in an analogous way for its computation and 
consider separately the two terms in ([42). We define the 



A and B term as the first (oc 



and the second 



(oc (1 — z)q ■ (q — zk)) terms of Eq. (|42|) respectively. 



4.3.1 A term 

The A termQ, owing to the k' 2 factor in the numerator, 
has no g || 1 collinear singularity. In addition, the z 2 factor 

2 This term has the same structure of the Cf term analyzed 
in Ref. Q, and can be simply recovered from the latter by 
replacing CpV gq — > C\P. The analysis that we perform here 
closely follows the one made in that paper. 



in the numerator causes a suppression of the integrand 
in the central region, and we can shift the lower limit of 
integration z cut in z down to zero: 



^ A = as ^ J dk h^(k) J dx f^(x)j\z z 2 P 



7r e q 2 (q - zk) 2 

We perform the rescaling q —. zl and use I as integration 
variable by substituting k' = k — zl, so that q — zk = 
z(l — k). Next we perform a simple fraction decomposition 
in order to separate the initial (i) state (g || 3 <==> 1 = 0) 
and final (/) state (1 || 3 ^==> I — k = 0) collinear 
singularities: 



l 2 (l-k) 2 l z + {l-k) 2 



I 



(l-k) 2 



(44) 



Beginning with the final state (/) collinear singular- 
ity, in terms of the new variables the A contribution to 
the jet cross section can be rewritten in the form 



dJ 



,l-2e 



dl 



7T e (Z - k)- 



da f / oft 
dJ 



i /.coll 

da^ 
dJ 



ckr /.finite 

dJ 



(45) 



which is particularly suitable for the analytic extraction 
of the divergences: the RHS contains three pieces, 1 ) the 
soft divergence, 2) the collinear divergence, and 3) a finite 
part. The integrand introduced in Eq. (Eq) is defined by 



7( 2 , i) 



^zP(z) [ dk fr(°>(fc) . MCa 

7T V 'J q V 'l 2 + (l-k) 2 



dx S { f } (k - zl, zl, xz; x)f^> (x) . (46) 



(0)/ 



The soft term in Eq. (^5|) is defined by evaluating the 
integrand in the soft limit z — > 0. In this lim it, t he jet 
distribution can be simplified by means of Eq. (29b), and 
leads to the constraint k 2 = E 2 . One obtains 



da^ soft 
dJ 



a. 



dz 

l-2e 



dl 



dk h^\k) 



10 z 

Ca 

S 

7T 

NC A 
l 2 + (l- k) 2 

C A 



7r e (i-fe) : 



rI(0,Z) 



dz 

y l-2e 



dl 



Ks(l - k) 



dxSf\k;x)f^(x) 



1 


7T 2 " 




2^2 


12 


m 



x dk dxhW(k)V< \k,x)f<?\x) , (47) 



where as usual, in the final result, we have collected some 
factors in such a way that they reproduce the LO jet ver- 
tex (0). The divergent factor exhibits single as well as 
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double poles, because our definition of the soft part in- 
cludes also the region where collinear and soft singularities 
merge. 

The pure collinear singularity can be isolated by eval- 
uating the integrand (^6|) in the collinear limit I = k, after 
having subtracted the soft term (I = fe, z = 0). The result- 
ing expression is clearly regular in the soft limit (z — » 0) 
and therefore contains a simple collinear pole. An UV cut- 
off A is introduced since the residue at the collinear limit is 



(2). 



no more integrable in the UV region. Thanks to Eq. (29c) 
the jet distribution simplifies to S 



da*/ 011 
dJ 



dl 



7T e (Z-fc) 2 



x [/(*, fe) - 1(0, fe)] 0(A 2 -(I- fe) 2 ) 



as 9A J dfe h (0) (k) 
x [zP{z) - 1] 



fe 2 
dl 



n £ (l-k) 2 
x J dxSf(k;x)f^(x) 

C A 



dz 

■A-2e 



0(A 2 -(I- fe) 2 ) 



[-- 




67" 


L 12e 


'S)' 


"36 



x / dx J dkh^(k)V^ 0) (k,x)f^(x) . (48) 

The remaining part is regular in the e — > limit and 
defines the finite term: 



^/.finite 

dJ 



1 dz 

z 



dl 



(l(z,k)~I(O,k))0(A 2 -(l-k) 2 )\ 
dz 



as £A f dk h (0) (k) f^^- 

* J q Jo (1-*) + 

x [(l-z)P(l-z)} f dxti°\x) 



dl 



MCa 



l 2 + (l- kf 



x Sf\l + z(k - I), (1 - z)(k - l),x(l ~ z);x) 



-V^(k,x)0(A 2 -I 2 ) 



(49) 



In the last equation we have performed the change of vari- 
able z — > 1 — z in order to simplify the expression with 
the introduction of the ()+ regularization for regularizing 
the 1/(1 — z) distribution at z = 1. 

Next we consider the term (i) with the initial state 
collinear singularity. We can write, in the same way as 
before, 



da] 
dJ 



,1—26 



dl 



I(z,l) 



_ da^ oft day 011 
~ dJ 



7Tj 2 
dJ + dJ 



j _ i, finite 

da X 



(50) 



where I is given by Eq. (|46|). It is trivial to see that the 
soft contribution, defined by evaluating the integrand I 
in Eq. (|5C|) at z — 0, is equal to the corresponding (/) 
term (14 71): 



da 1 



da f / oit 



(51) 



dJ dJ 

As to the collinear piece, we note that in the collinear lim it 
1 = the jet distribution reduces (by applying Eq. (29e)) 



(2) 

to Sj , and one gets, using always the same UV cutoff A, 
performing the I integration and changing the integration 
variable z — > 1 — z, 



da_^ 
dJ 



i, coll 



dz 



dl 



^ / ttJ 2 



x [/(z,0) -1(0,0)] 6>(yl 2 -I 2 ) 
a 9A J dfe h (0) {k) J dxf (0) {x) 

f 1 fA 2 Y f 1 dz 



v £ V^V Jo (1 - z )+ 
x U(°)(fc,a;z) + 2 



(l-z)P(l-z) 



/ dz 


"ln(l - z)" 


[ (l-«) J 


Jo 



x (1 - z)P(l - z)K (0) (fc 



(52) 



The left contribution is regular in 4 dimensions and 
defines another finite term 



^finite 

dJ 



°»/f / Jrl/CM)-.r(o,Q 

-(/(z,O)-/(O,O))0(^ 2 -i 2 )] 

dz 



= a^jdH h ^k)f o - 



z). 



x [(l-z)P(l-z)} dx6°\x) 



dl 



l z + (l-k) 2 



x 5 ( 7 3) (fc - (1 - z)Z, (1 - z)Z, x(l - z);x) 



V^ o) (k,xz)0(A 2 -I 2 ) 



(53) 



4.3.2 B term 

The last piece to be analyzed is the B term, which reads 
% B = as ^jdkh^{k)jdxf^{x) 

x f dz (l-z)P(z) ' <lk 



7T F k' 2 



x q 2( iq -tl sf\k', q ,xz; X ) 
q z (q — zky ' 

da^ 11 derLL da c £ nst 

dJ + dJ + dJ 



(54) 
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and will be decomposed in a collinear divergent piece, an 
energy dependent piece which is just the LL contribution 
to the cross section coming from gluon emission in the 
central region, plus a finite and constant in energy term 
as shown in the last line of the above equation. It can be 
noticed immediately the presence of a g || 1 collinear sin- 
gularity corresponding to the k = pole. As already com- 
mented for the similar case there are no other singular- 
ities in the fc'-integration (neither q = nor q — zk = 0) 
except for z — » corresponding to gluon 3 being in the 
central region, which gives the LL contribution. Here we 
expect a soft singularity needed to cancel the one in the 
gluon trajectory ([l6|). 



easily obtains 



We recall that the jet distribution functions become 
essential in disentangling the collinear singularities, the 
soft singularities, and the leading logs pieces. 

The basic mechanism are the same as for the quark 
case (remember here that gluon 3 is playing a special role 
because of the manipul ation in the integrand performed 
at the beginning of Sec. 4.3): 



When the outgoing gluon 1 is in the collinear region of 
the incoming gluon g, i.e., y\ — > oo, it cannot enter the 
jet; only gluon 3 can thus be the jet, y 3 is fixed and 
no logarithm of the energy can arise due to the lack 
of evolution in the gluon rapidity. No other singular 
configuration is found when J = {3}. 
In the composite jet configuration, i.e., J = {1,3}, the 
gluon rapidity is bounded within a small range of val- 
ues, and also in this case no log s can arise. There could 
be a singularity for vanishing gluon 3 momentum: even 
if the 1 || 3 collinear singularity is absent, we have seen 
that, at very low z, a soft singular integrand arises. 
However, the divergence is prevented by the jet cone 
boundary, which causes a shrinkage of the domain of 
integration ~ z 2 for z — > and thus compensates the 
growth of the integrand. 

The jet configuration J = {1} allows gluon 3 to span 
the whole phase space, apart, of course, from the jet 
region itself. The LL term arises from gluon configu- 
rations in the central region. Therefore, it is crucial to 
understand to what extent the differential cross sec- 
tion provides a leading contribution. It turns out that 
the coherence of QCD radiation suppresses the emis- 
sion probability for gluon 3 rapidity j/3 being larger 
than the rapidity y\ of the gluon 1, namely an angular 
ordering prescription holds. This will provide the final 
form of the leading term, i.e., the appropriate scale of 
the energy and, as a consequence, a finite and definite 
expression for the one-loop jet vertex correction. 



Let us isolate in Eq. (Q) the initial state g || 1 collinear 
singular contribution and define, as usual, the collinear 
term by setting k' = (except in the 1/k' 2 pole), and by 
introducing an UV cutoff. Observing that the jet distri- 

1— 1 tn\ 

bution, because of Eq. (|29d[ ), reduces to Sj (k;xz) one 



/i 0) (*) 



xz) 



(55) 



The LL part can be extracted exactly with the same 
procedure already followed in Q, thanks to the suppres- 
sion of the partonic cross section when gluon 3 is emit- 
ted at larger rapidity w.r.t. that of gluon 1. Infact, when 
gluon 3 is in the central region, gluon 1 must be the jet. In 
this case, the azimuthal average of the cross section w.r.t. 
gluon 3 azimuthal angle 03 at fixed gluon 1 momentum 
yields 



(1 - z)q ■ (q - zk) 
q 2 (q — zk) 2 



\o(E 3 - z{Ex + E 3 )) 



(56) 



This relation is exact and clearly shows that, outside the 
angular ordered region 



E 3 > E l 
z 1 — x 



h > 



2/3 < Vx , (57) 



there is no contribution to the cross section. In practice, 
by taking into account the variation of h^\k) during the 
averaging procedure, instead of a strictly vanishing con- 
tribution we have a strong suppression. Note that, in the 
limit q — ► (which includes the soft region], the variation 
of k' goes to zero as well, so that Eq. ([57|) is really an 
accurate statement in the "dangerous" part of the phase 
space. Moreover, Eq. (5(j) shows that the 1/q 2 kinematic 
dependence of the LL kernel governs the differential cross 
section up to the very end of the angular boundary. 

Therefore, we define the LL contribution in the "lower 
half region" y' 3 > by 

:= a s [ dk h^(k) f dk' ^ J^(°)(fc') 
dJ J J it Tr £ q 2 8 

f ^O^-ziE.+E,)) 
x f dxSf\k',x)f^(x) 



(58) 



d.T J dk J dk' h^{k)K^^\k,k') 



x log 



Ej + E 3 g 



vW(k',x)fQ\x), 



having imposed the jet condition J = {1}. 

The remaining part of the B term is finite in 4 dimen- 
sions and constant in energy, so that we can set e = and 
^cut = to define the constant part 



dfJ const 

dJ 



da B 
dJ 



da c ° 11 
dJ 



da^ 
dJ 



2cut = • 

s = 



(59) 



10 



J. Bartels et al.: The NLO Jet Vertex for Mueller-Navelet and Forward Jets: the Ghion Part 



5 The NLO jet vertex: sum of real and virtual 
corrections 

Having completed the calculation of both the virtual and 
real corrections in the whole phase space, we are going to 
collect all partial results and to show that the complete 
one-loop jet cross section can naturally be fitted to the 
form of Eq. @. Table [3] summarizes the decomposition 
of the one-loop jet cross section and gives the references 
of the various contributions. 

In Sec. |§| we have presented the virtual contributions 
to the jet cross section which, after renormalization of the 
coupling, assume the form of Eq. (fi"5|). The contribution 
coming from the 27 q impact factor correction has been 
combined with the "upper half region" real contribution 
to give the full impact factor of the upper quark q . 

The remaining virtual terms and the LL contributions 
of both the negative and positive rapidity regions (the lat- 
ter contribution given in ( f58| ) ) can be conveniently rewrit- 
ten in the form 



dcr ( virt ) 
dJ 



dcr ( virt ^ 



+ dj 



da( r 
dJ 



LL 



a s J dx J dfcdfc' h^(k)K w (k,k') 



x log 



s (k,k ) 



(0)/ 



(60) 



Ca 
n 



k 1 



12 



C A 1 IIC A - 2N f 1 



7r e z 
67\ 5 Nf 



()7T 



•S(i 



18 7T 



k 

b log -a 
M 



x dx dk h^(k)V^(k,x)f^(x) , (61) 



where the kernel K^°> is the full BFKL kernel and the 
energy scale sq in Eq. ( |60| ) has been chosen as 

s (k, k') := (|fe'| + | 9 |)(|fe| + \q\) = (Ej + E 3 )(E 2 + E 3 ) . 

(62) 

Let us stress that the energy scale in Eq. ([32]) arises 
naturally when one requires impact factors and PDFs to 
have standard collinear properties and the remaining non- 
leading-log term to be finite in both the physical e — > 
and high-energy s — > oo limits. Choosing a different 
scale of the energy requires the introduction of additional 
NLL operators, which has to be added as multiplicative 
corrections to the Green's function. If, for i nsta nce, we 
adopt so = |fc||fc'|> then the Green's function ( |l0b| ) has to 
be replaced by 

G(xs,k,k') = (l + a,H L ) l + a s i^°>log-^ 

x (1 + a s H R ) (63) 
H L (k,k') = -K^(k,k')log l -^±l^ 
H R (k,k')=H L (k',k) . (64) 



We now consider the sum of all e-divergent contribu- 
tions coming from both real and virtual corrections. All 
double poles coming from (pl h , f [47|) and f[5l"l) cancel out. 
The single poles come from~^|), @, (|§, ©,(11), 
(p5|). They combine to give exactly the DGLAP splitting 
functions terms which define the PDF one-loop correc- 
tions. One has only to note that 



(l-z)P(l-z) 
(1-*)+ 



P(z) 



(1-*)h 



+ z(l-z) 
(65) 



and that the P g , 

by 



LO DGLAP splitting function is defined 



P gg (z) - 2C A 



1 



z 

HC A 



(i- z y 



+ 2(1 - Z) 



(66) 



The non-LL total singular contribution adds up to 



da; J dk h^\k) 



da singular a s \ / A 2 V 

dJ ~^e\Jj?) 

+vW(k,x)[P gg ®fW](x)} . (67) 



The cutoff A introduced in all the collinear subtraction 
has the physical interpretation of factorization scale and 
can therefore be identified with (if. Eq. (|67j ) shows that 
the only singular terms remaining after the sum of real 
and virtual corrections are of collinear origin, and actually 
amount to the expected collinear singularities stemming 
from parton radiation out of incoming gluons. They can 
be absorbed in the PDFs according to Eq. ( |l3| ) ; of course 
only the term b = g in the sum is reproduced. 

Finally all the finite and not LL terms from (J6l| 

©, dp, ©, ©. ©> ©, ©> ©, ©. © 

and (|59|) contribute to the jet vertex 




d(T finite 

dJ 



a s dx dkh^(k)V g ^(k,x)f^(x) , (68) 
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Table 1. Schematics of the decomposition of real and virtual one-loop corrections to gq scattering and references of the 
corresponding equations. 



virtual 


real qgg 


real qqq 


oj w n q n q 


Vz <0 


J/3>0 






A 


73 


Cf 




/ \ 

soft 

LL coll h ( q 1} LL 
(|SG|) del]) (26b) (26a) 


f 


i 


coll LL const 

(R) (R) 


coll fin 


coll fin 

d39|) (Fob 


soft coll fin 
d«t ® (pj) 


soft coll fin 
P P (p|) 



which defines the NLO correction to the gluon-initiated 
jet vertex 

6 IT 3 7T / vl 2 4 f 36 7T 



fe 2 

-b log -J 
A 1 



K (0) (fc, 



+ / dz ^(fe.a*) [^§^(1 -z) 

J L 7T 

+ 2^(1 - z P 1 - z) ^ 

7T \ 1 — Z 



TV/ /" dfc' 



dzP qg (z) 



2 + fc 72 5 J (*,g,a;«;s) 



--^6)(^-fc' 2 )T/ q (°)(fc,xz) 
\o(A 2 -{q-zk) 2 )sf \k,x) 



Ca 


Jo ( 


IT 






dl 


'J 









dz 



{(l-z)P(l-z)} 



ttP \l 2 + {l-k) 2 
S { j\zk + (1 - z)l, (1 - z)(k - l),x(l - z);x) 
+ Sf } (k - (1 - z)l, (1 - z)Z,.t(1 - z);x) 

- 0(A 2 - I 2 ) [Vg (0) (fc, a:) + (k, a?z)] | 

q • (g - zfc) 



/"dfe 



P(z) (1-z) 



q 2 (q — zk) 2 



x /i(°)(fc')^ 3) (fc',q,a;z;a;) + 
- ^0{A 2 -k' 2 )V^{k,xz) 

-^0{\q\ - z(\ q \ + \k'\))V^{k\x) 
zq 



(69) 



Like for the quark case the gluon originated vertex clearly 
depends on the jet definition and on three scales: the en- 
ergy scale so (via the subtraction of the LL term oc 1/z), 
the factorization scale A — fip and the renormalization 
scale /i. 

6 A consistency check: jet vertex versus 
impact factor 

In this section we compare the structures of the high en- 
ergy factorization formulae in the two cases of jet cross 
section — discussed in this paper and in [|] — and of in- 
clusive partonic cross section — presented in Jl3],[l4) . This 
is interesting both for understanding the different struc- 
ture of the singularities involved in the two formulations 
and, on the other hand, to underline the equivalence in 
the LL subtraction procedure which will allow us to use 
the results of two-loop calculations for the generalization 
to higher orders presented in Sec. 0. 



6.1 Jet cross section versus inclusive cross section 

Let us consider the quark-hadron one-jet inclusive factor- 
ization formula 



da 
dJ 



qH 



J2 J dxdkdk' h q (k)G(xs,k,k')V 3 (k',x)f a {x) 

(70) 



(cfr. Eq. (10a)) and compare it to the one describing the 
inclusive quark-parton cross section 



da™ cl = / dfcdfc' h q (k)G(s,k,k')h a (k') . (71) 



Formally, one can obtain the former from the latter by 
simply replacing the partonic impact factor of the lower 
incoming parton a with the product of jet vertex and PDF 
— we refer to this product as "jet impact factor" . There is 
nevertheless an important difference between the one-loop 
correction to the impact factor and the one to the "jet 
impact factor" [V * = * / W + * f {1 \ In the 
gluon-initiated process, the jet impact factor contains the 
full collinear singularities (included in f^ 1 ') relative to the 
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g — > q and g — > g splittings, as shown in Eq. (|67|). The 
gluon impact factor, on the other hand, contains the full 
g — * q collincar singularity but only a part of the g — > g 
one. In fact, the residue of the 1/e pole of the gluon impact 

factor correction is J dz [P qg (z) + P gg (z) — 2Ca/z}. 
This is because the subtracted term 2Ca/z has been taken 
into account in the LL contribution (see Eq. ([58])). The 
same is true for quark-initiated processes, the residue of 
the 1/e pole in the quark impact factor correction hq 
being Jdz [P qq (z) +P gq (z) - 2C F /z}. 

In the jet case, it has been possible to factorize the 
full collincar singularities into the PDF because of the jet 
distribution, which selects only some final states out of 
the whole phase space. 



(1,0) 




(i-z, *') 



Fig. 5. LAB frame view of the quark + reggeon —* quark 
+ gluon blob involved in the quark-initiated jet cross section 
(in particular the Feynman diagram responsible for the 1 || a 
collincar singularity yielding 1/e x P gq is shown) . 6i and 83 de- 
notes the outgoing quark and gluon polar angles respectively. 
In parentheses close to each particles are the longitudinal mo- 
mentum fraction (w.r.t. the incoming parton) and the trans- 
verse momentum. 

How does this happen? Let us discuss this issue in more 
detail for the simpler case of quark-initiated scattering. 
Fig. ^ shows a possible configuration of the quark-reggeon 
blob qg* — > qg. In the inclusive cross section, all config- 
urations contribute to the cross section. When #3 < 9%, 
the contribution is entirely assigned to the impact factor 
correction. However, when #3 > 9±, according to angular 
ordering, the part oc 2Cp/z of the squared matrix element 
is interpreted as LL term. Since the collinear region which 
should provide the P gq splitting function times the 1/e 
pole is just 9i ~ 0, it is clear that it is always affected by 
the LL subtraction, and therefore the impact factor correc- 
tion contains only the "non singular" part P gq (z) — 2Cf/z 
of the splitting function. 

In the jet case, in the collinear region 9\ ~ 0, we per- 
form no LL subtraction, because the gluon is constrained 
to be the jet, hence the 1/e collinear pole multiplies the 
full P gq splitting function (and the LO jet vertex). The LL 
subtraction is performed when 9 3 > 9i and the jet con- 
sists of only the quark. In this case there is no collinear 
singularity associated with the LL subtraction. 

6.2 Relation between M 1 ) and V" (1) 

One could think to check the consistency of our calcu- 
lation by integrating the factorization formula of the jet 



observables with respect to the jet variables in order to 
recover the quark-hadron inclusive cross section 

d<J qH = J2 [ dxdkdk> h q (k)G(xs,k,k')h 3 (k')f^(x) , 

(72) 

obtained by convoluting the quark-parton inclusive cross 
section ( f7l| ) with the bare-] parton densities (x) . 

This is not the case, because the single-jet configura- 
tions would be counted twice. Consider, for instance, the 
configuration of Fig. |L It would be counted once when the 

quark is the jet and a second time when the gluon is the 

(3) 

jet. On the contrary, if we replace the jet distribution Sj 

— (3) 

with another distribution S 'y with the property that 

J dJ S { f\p 1 ,p 3 ,^;x) = 1 , dJ:= dyjdPjd^j , 

(73) 

(i.e., a function that counts each final state once and only 
once), and define a "modified jet vertex" obtained 
from l/W by replacing Sj 3 ^ — > sf \ then the integration 
over the jet variables of the "modified one-jet inclusive 
cross section" yields the inclusive cross section ([72]), be- 
cause 

J dJ V * f = h* . (74) 

More precisely: 

x j dz [V a ( - 0) (k,xz;J)P aq (z) + V a {0) {k,xz;J)P ag {z)]^ 
= hP(k), (a = q,g). (75) 

Despite of the difference between the partonic and jet 
impact factors, the remaining factors in the respective 
factorization formulae are the same. In particular, the 
Green's function at scale |fe||fe'| in the form of Eq. ( |65| ) 
is identical to the one defined in |l3| for the partonic case, 
i.e., it shares the same impact kernel Hl and Hr. 

This is due to the method adopted for the definition of 
the partonic impact factor (resp., the jet vertex), i.e., for 
the separation of the central region from the fragmenta- 
tion region. In practice, whatever the selection of the final 
state for the observable under study, the precise contri- 
bution of the fragmentation region of each of the incom- 
ing particles has been defined according to the following 
scheme: 

- take the expression of the squared matrix element (sup- 
plied eventually with the distribution S providing the final 
state selection) valid when the gluon is in the fragmenta- 
tion region of the incoming particle; 

3 One must use the bare PDFs because the collinear singu- 
larities are embodied in the partonic impact factors. 
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- extract the central region limit from the above expres- 
sion. It should be proportional to the BFKL kernel times 
the LO impact factor times an enhancing factor describing 
the large extension of the central region (in our variables, 
this enhancing factor is 1/z), provided the final state se- 
lection is inclusive in the central region. This is the case 
for the one- and two-jet inclusive observables that we are 
studying; 

- at fixed quark momentum, subtract from the fragmenta- 
tion region expression the one obtained as central region 
limit. This subtraction has to be performed according to 
the angular ordering prescription #3 > 9\ 1/3 > yi 
in the soft gluon region. 

This defines on one side the impact factors, and on the 
other side what is not impact factor, i.e., what contributes 
to the Green's function. The universality of the Green's 
function (at least at one-loop) should now become clearer: 
it is due to the same scheme of defining the LL subtraction, 
yielding in all cases the BFKL kernel times the LO impact 
factor. The details of the angular ordering out of the soft 
gluon region and the choice of the energy scale determines 
the form of the impact kernel Hl.r- 



7 Factorization formula at higher orders 

In this section we discuss some issues related to the gener- 
alization of the factorization formula (|l0| ) to higher order 
perturbation theory. 

It was the main purpose of this paper (and of the pre- 
vious one [E]) to provide a factorization formula for the jet 
production cross section at high energies, accurate in the 
next-to-leading log s approximation. To this end we have 
studied the jet cross section formula in next-to-leading or- 
der, 0(af). The leading logarithmic approximation (LL) 
retains the coefficient of Ins. What we have computed in 
this paper is the next-to-leading logarithmic approxima- 
tion (NLL), i.e. the terms without the Ins enhancement. 
The generalization to higher order in a s therefore con- 
tains, in LL, terms of the form a™ (In s)™ -2 , and in NLL 
contributions of the form a™(lns)™ -3 . We are interested 
in summing both the LL and NLL terms to all orders a s . 

In generalizing our result to all orders we have to rely 
on Regge factorization properties of scattering amplitudes 
in the high energy limit. Regge theory states that ampli- 
tudes (both for elastic scattering, and also for the inelastic 
production of particles in the multiregge limit, when de- 
scribed by the exchange of Regge poles) obey factorization 
properties. In QCD the gluon has been found to reggeize, 
and, as it is discussed in some detail in the original BFKL 
papers jB and in the NLO calculations of the BFKL ker- 
nel [^scattering amplitudes with color octet exchange 
are expected to satisfy these factorization properties. In 
QCD calculations, this leads to the 'Regge ansatz' for elas- 
tic parton-parton scattering and for inelastic gluon pro- 
duction in the the multiregge limit. Strictly speaking, this 
ansatz is a hypothesis, suggested by Regge theory and by 
the observation that the gluon reggeizes. The strongest 
support for the correctness of the assumption comes from 
the result that the LL BFKL equation, derived from the 



Regge ansatz, in the color octet channel satisfies the boot- 
strap condition. More recently, the Regge ansatz has been 
used also in NLO calculations the bootstrap con- 

dition in the NLO BFKL equation was first formulated 
in its weak form inside elastic scattering amplitudes and 
after in its strong form, which is the one necessary for the 
self-consistency of the assumption of Reggeized form of 
the production amplitudes |l^,|l8|. Very recently the va- 
lidity of the bootstrap condition (strong form) in the NLO 



BFKL equation has been proven in 1 19 

In our case, we have used the Regge ansatz for the 
factorizing form in (|Io|): in the jet production amplitude 
the large rapidity gap between the outgoing quark (at 
the upper end) and the jet (at the bottom) leads to the 
exchange of a reggeized gluon, and the upper and lower 
vertex become in dep endent from each other. In the cross 
section formula (|10|) this factorization property leads to 
the impact factor and the jet vertex. In our calculation, 
it has been one of our tasks to show that the jet vertex 
fits into this factorization scheme: when integrated over 
the jet variables, the square of the "modified jet vertex" 
turns into the gluon (or quark) impact factor which has 
been studied before. Without the integration over the jet 
variables we have, in principle, a new situation: we find 
that Regge factorization still works at the NLO level. One 
expects that this factorization continues to hold also in 
higher order a s , i.e. the impact factor at the upper end of 
the gluon and the jet vertex allow for a systematic expan- 
sion in powers of a s . Similarly, the exchanged gluon which 
in lowest order is elementary turns into a reggeized gluon 
with a trajectory function being given as a power series in 

When generalizing to higher orders, i.e. to the produc- 
tion of more gluons, we follow the strategy pursued in the 
NLO calculations of the BFKL kernel: the Regge ansatz 
for multi-gluon production amplitudes states that vertices 
that are separated by large rapidity gaps are independent 
of each other, i.e. they factorize. Before turning to the all- 
order formula, it may be instructive to take a closer look 
into the next order correction to the cross section formula, 
the two-loop jet cross section 0(a^). 

Since final states may consist of four, three, or two par- 
tons we have terms of order log 2 s and a* log s (terms 
without Ins belong to NNLL and will not be considered). 
Regge factorization leads us to expect: 

i_^ (2) = I log 2 ff h (0) K (0) K (0)y(0) /( 0) 

af dJ 2 sr 

+ l og ™\ h W K (v) V (°) f (0) 

SR L 

+ (H L K<V + K^H R + kW)vW / (0) 
+ h^K^V^ / (0) + h (0) £r(°V (0) 

(76) 

The first term on the RHS is the LL term given by the 
iteration of the leading kernel K^°\ The square brackets 
collect the NLL contributions: the impact factor correc- 
tion, the contribution from the products of impact kernel 
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Fig. 6. Schematics of amplitude diagrams involved in the calculation of the partonic cross section at relative O(a^): (a) two 
particle final state with two-loop virtual corrections; (b) three particle final state with one-loop virtual corrections; (c) four 
particle final state at tree level. The final state particles are ordered with increasing rapidity from the top to the bottom. 



and leading kernel, the NLL kernel, and finally the correc- 
tions to the jet vertex and to the PDF. 

In the second step we relate the terms of Eq. ( |76|) to 
the various final states and kinematic regions contributing 
to the 0(af) cross section. We have three kinds of final 
states: two-parton final states (with two-loop virtual cor- 
rections), three-parton final states (with one- loop virtual 
corrections) and four-parton final states at tree level. Be- 
cause of the reggeization of the gluon, there is no one-to- 
one correspondence between the number of rapidity gaps 
and the power of Ins. Let us go through the final states 
in some detail. 

Consider first the four-parton final state of Fig. |(]c. It 
can be shown that, when the two central gluons are emit- 
ted at large subenergies S> tj, i.e., final state particles 
are separated by large rapidity gaps, the cross section is 
described by the iteration of two BFKL kernels (actually 
their real parts), one for each emitted gluon. Since the 
external partons are isolated in rapidity space, they are 
described by LO partonic impact factors in the inclusive 
process, and by the LO jet vertices in the jet process. 

NLL contributions arise from quasi-multi-regge kine- 
matics, where one (and only one) pair of particles is emit- 
ted at small subenergy (small rapidity gap). If the pair 
with small subenergy is formed by the two central gluons 
3 and 4, there is no ambiguity to associate the correspond- 
ing contribution to the NLL BFKL kernel This is 
one of the so called "irreducible contribution" of Ref. || , 
because it can directly be interpreted as contribution to 
regardless of the definition of the impact factors. 
Note that the rapidity intervals y\ — 2/3 and 2/4 — 1/2 must 
be large in order to ensure the contribution to be NLL. 
Therefore the external particles 1 and 2 are still described 
by LO impact factors or jet vertices. 

When one of the central gluon approaches an external 
particle, say gluon 3 enters the fragmentation region of 
particle 1, S13 and 2/1 — 2/3 become small. In this case we 
have the contribution to the jet vertex and PDF correc- 
tions. The definition given in the one-loop situation still 
works here, always because of high energy factorization: 
since S34 and S42 are much larger than S13, the emission 
of gluon 4 is described by the real part of the LL BFKL 
kernel, independently of what is going on in the upper 
or lower vertices. Therefore, we can apply the machinery 
of the one-loop calculation factoring out the b — > 2 + 4 



subprocess (described by h^K^, and consider only the 
blob a + g*(fc) — > 1 + 3 as done before. What one obtains 

— after the subtraction of the central region contribution 
of gluon 3 according to the angular ordering prescription 

— is the sum of the one-loop (real) corrections to the ver- 
tex and to the PDF, the impact kernel Hr, and a residual 
contribution that can only be assigned to the NLL kernel 

This is one of the so called "reducible" contribution, 
in that it depends on the definition of the "(jet) impact 
factors" . The key point is that also the reducible contri- 
butions to K^ 1 ' are the same as those obtained in the 
inclusive case, because they are related (just like Hr) to 
the subtraction of the central region from the fragmen- 
tation region, and we have shown in Sec. |^ that the LL 
subtraction is defined exactly in the same way in the jet 
cross section as well as in the inclusive cross section. 

The treatment of the phase space region where gluon 
4 is close to parton 2 is equivalent to that adopted in the 
inclusive process and provides the upper impact factor 
correction the other impact kernel Hl and another 
reducible contribution to K^ 1 '. 

A similar analysis can be performed for the three par- 
ton final state emission, Fig. ^b. When the gluon is in 
the central region, one has the product of the real part 
^rcai °f the- kernel (gluon emission) and of the virtual part 
(provided by the leading virtual corrections), but 
also virtual corrections to the partonic or jet impact fac- 
tors (given by virtual corrections attaching only to the cor- 
responding partonic line) and an irreducible contribution 
to (given by virtual corrections attaching to the out- 
going central gluon and to the exchanged gluons). When 
the emitted gluon approaches the lower partonic line, the 
contribution to the cross section can be at most NLL, and 
this occurs only if there are virtual corrections "filling" 
the large rapidity gap 2/3 — 2/2 providing a factorized one- 
loop gluon trajectory (virtual part of the BFKL kernel). 
Therefore, also in this case we can factor out the upper 

blob b — > 2 (described by h^K^t an d apply the proce- 
dure of separation of central region and fragmentation re- 
gion to the lower blob a + g*(fc) — > 1 + 3 as done before. We 
obtain again — after the subtraction of the central region 
contribution of gluon 3 according to the angular ordering 
prescription — the sum of the one-loop (real) corrections 
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to the vertex and to the PDF, the impact kernel Hr, and 
another reducible contribution to 

The emitted gluon approaching the upper partonic line 
is treated in the same way of the inclusive case, giving 
again a correction to the upper partonic impact factor, the 
composition of the impact kernel Hi with the virtual part 
of the BFKL kernel and yet another reducible contribution 
to A«. 

The two-parton final state contribution (see Fig. ^|a) 
is the simplest to analyze, since the jet has a trivial struc- 
ture, namely it consists only of the lower outgoing quark, 
and is described by the LO jet vertex oc h^°\ The 
comparison with the partonic case is therefore straight- 
forward and we obtain: the iteration of two gluon trajec- 
tories a/ 1 ) (to be included in K^K^) 7 the product of 
the virtual correction to the impact factor and the gluon 
trajectory (to be included in hP^K^), the product of the 
trajectory and the virtual corrections to the lower quark- 
reggeon vertex (to be included in K^V^ and f^ 1 ') 
and the two-loop gluon trajectory uj^ (the last irreducible 
contribution to K^- 1 '). 

From Eq. ( |76| ) and from the subsequent discussion one 
sees that the quantity governing the high energy behaviour 
of the higher order processes at NLL level is the BFKL 
kernel /C, which has the perturbative expansion 



(77) 



Its NLL term has been determined for inclusive par- 
tonic scattering. In a perturbative expansion of the cross 
section, the NLL kernel AT' 1 ' enters the first time at rela- 
tive order a^, i.e., at two- loop. 

The most straightforward way of generalizing to all or- 
ders in a s is to adapt the strategy used in the NLO study 
of parton parton scattering to our case of one-jet-inclusive 
processes. The main difference between the two cases is 
that in the latter the impact factor of the lower incoming 
parton is replaced by the convolution of the PDF with 
the jet vertex, as was already shown in the one-loop cal- 
culation. Therefore, we keep the general structure given in 
Eq. (10a), together with the perturbative expansions ( [TT| ) 
of the energy independent factors. Only the Green's func- 
tion G needs to be extended to NLL level. 

Before we write down the general all-order formula, let 
us return, once more, to the question of the energy scale 
so- For future applications it will be convenient to use the 
symmetric Regge-type energy scale sr:= |fci||fc 2 |. As we 
have discussed before, our requirement of giving the im- 
pact factor and the PDF's the correct collinear properties 
has led to the energy scale (62), and a change of the 
energy scale will be induced by the additional operators 
(1 + a s Hft), (1 + a s Hi). As discussed in (]]J, this leads 
to the following form of the NLO BFKL Green's function 
for the energy scale sr := |fci||fc 2 |: 



G(xs, fei,fe 2 ) = J 



duj f xs' 
27ri \sr / 

x {k x \{\ + a s H L )[u - K]-\l + a s H R )\k 2 ) . (78) 



in terms of the impact kernels (|64|) and of the BFKL ker- 
nel (|77j). Inserting this Green's function into (|l0a[) and 
expanding in powers of a s one readily reproduces J7q) . 

We conclude with a final formula for the Mueller-Na- 
velet jets: 



^ 2 
dJ 



— - = y~W da;idx 2 / dfcidfc 2 f a (xi)V a (ki,xx) 

' " a,b 

G(xix 2 s, fei, fe 2 )Fb(fe 2 , x 2 )f b (x 2 ) 



(79) 

where a, b = q, g, and the subscripts 1 and 2 refer to jet 1 
and 2 in hadron 1 and 2, resp. All elements are known in 
LO and in NLO; in particular, the jet vertex expressions 
), and Vq (1) from (105) of §. 



vS from 



8 Concluding remarks 

In this paper we have completed the analytic part of the 
NLO corrections of the jet vertex which appears in the 
cross section formulae for Mueller-Navelet jets at hadron 
hadron colliders and for forward jets in deep inelastic elec- 
tron proton scattering. The final result of our study is 
summarized in eq. (|79|), the analytic expression of the 
cross section for the production of Mueller-Navelet jets 
at the Tevatron or at the LHC. 

Apart from the interest in performing a consistent NLO 
analysis of the BFKL Pomeron at hadron hadron colliders 
or in deep inelastic ep scattering, the results of our anal- 
ysis may rise some theoretical interest. Due to the very 
special kinematics of the jet production processes, the jet 
vertex lies at the interface between two different high en- 
ergy limits, the hard scattering regime and the Regge limit 
( small- X limit): a priori it was not clear whether, at NLO 
accuracy, it would be possible to separate the collinear di- 
vergences from the BFKL-type gluon production. We have 
obtained an affirmative answer. Previous experience with 
partonic impact factors has served as a valuable guide in 
performing this separation inside the jet vertex. 

What remains is the numerical evaluation of the cross 
section formulae, using the results derived in this paper. 

As the first step, one has to specify the jet functions Sj\ 
i.e. we have to decide on a jet algorithm. In addition, ex- 
perimental cuts have to be formulated. We hope to be able 
to report first numerical results in the near future. 
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